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ZEN2: A narrow J-band search for z ~ 9 Lya emitting 
galaxies directed towards three lensing clusters. 
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ABSTRACT 

We present the results of a continuing survey to detect Lya emitting galaxies at red- 
shifts z ~ 9 : the ZEN ("z equals nine") survey. We have obtained deep VLT/ISAAC 
observations in the narrow J-band filter NB119 directed towards three massive lens- 
ing clusters: Abell clusters 1689, 1835, and 114. The foreground clusters provide a 
magnified view of the distant universe and permit a sensitive test for the presence 
of very high-redshift galaxies. We search for z ~ 9 Lya emitting galaxies displaying 
a significant narrow-band excess relative to accompanying J-band observations that 
remain undetected in HST/ACS optical images of each field. No sources consistent 
with this criterion are detected above the unlensed 90% point-source flux limit of the 
narrow-band image, Jnb = 3.7 x 10~ 18 ergs~ 1 cm -2 . To date, the total coverage of 
the ZEN survey has sampled a volume at z ~ 9 of approximately 1700 co-moving 
Mpc 3 to a Lya emission luminosity of 10 43 ergs~ 1 . We conclude by considering the 
prospects for detecting z ~ 9 Lya emitting galaxies in light of both observed galaxy 
properties at z < 7 and simulated populations at z > 7. 
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1 THE SEARCH FOR THE MOST DISTANT 
GALAXIES 

Observations of distant galaxies provide a direct view of 
the early stages of galaxy evolution as well as probing 
the physical conditions of the high-redshift intergalactic 
medium (IGM). The advent of the Advanced Camera for 
Surveys (ACS) deployed on the Hubble Space Telescope 
(HST) in addition to wide field optical cameras oper- 
ated on 8m class telescopes has provided access to rela- 
tively large samples of galaxies located at z ~ 6 and be- 
yond. Observations of distant Lyman drop-out galaxies (e.g . 
IStanwav. Bunker fc McMahonl 120031 ; iDickinson et all 120041 ') 
and Lya emitting galaxies selected via narr ow-band pho- 
tometry (|Rhoads et al.ll2003l ; IHu et al.ll2004h indicate that 
these bright, distant galaxies cannot be the sole agents of 
the global re-ionisation demonstrated by stu dies of Gunn- 
Peter son absorption in high-redshift quasars l|Bunker et al.l 
2006). Current solutions to this dilemma centre upon the 
potential contribution of faint galaxies at z ~ 6 - pointed to 
in the very deepest, yet numerically smallest, ACS samples 
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- or the possibility that an earlier epoch of more intense star 
formation was responsible for the observed re-ionisation. 

Perhaps the most important contribution to the debate 
concerning the physical nature of bright, high-redshift galax- 
ies has come from Spitzer space telescope obs ervations of 
rest-frame optical em ission in these systems fe.g. lEvles et al.l 
l2005l : lYan et alj|2006h . Optical to near infrared (NIR) obser- 
vations of z ~ 6 galaxies sample rest-frame emission blue- 
ward of the 4000A discontinuity and are thus sensitive pri- 
marily to the spectral contribution from younger stellar pop- 
ulations. The addition of Spitzer/IRAC bands, particularly 
at 3.6|r.m and 4.5u.m, samples the potential contribution of 
older stellar populations in these galaxies. From the limited 
number of bright z ~ 6 galaxies studied to date, there has 
emerged a picture of their stellar populations as being rela- 
tively old (up to 700 Myr) and massive (up to 3 x 10 10 M Q ). 
The extrapolation of these integrated star formation histo- 
ries to earlier cosmic times points to an epoch of potentially 
intense star formation in the predecessors of bright z ~ 6 
galaxies extending to redshifts z ~ 10. 

Though compelling evidence points to the existence of 
actively star forming galaxies at z > 7, the direct obser- 
vation of such sources is far from straightforward - mainly 
due to the extreme faintness of high-redshift galaxies. The 
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brightest galaxies observed at z ~ 6 display total AB 
magn itudes of order 24 (e.g. IStanwav. Bunker fc McMahonl 
2003; Hu et al. 2004). Galaxies at z > 7 - including cur- 
rent samples of candidate systems - can be reasonably ex- 
pected to display sig nal levels AB > 25 in NIR wavebands 
(Bou wens et alf2004h . Obtaining a spectroscopic redshift for 
such faint, continuum selected systems with currently avail- 
able technology is challenging (though not impossible, c.f. 
iKneib et alj l2004h . Studies employing narrow-band filters 
or long-slit spectral observations are sensitive to the sub- 
set of high-redshift galaxies that emit a significant fraction 
of their energy in the form of narrow spectral lines - typi- 
cally the Lya emission line. Deep narrow-band imaging and 
subsequent spectroscopic observations have been employed 
successfully to generate high spectral c ompleteness sam ples 
of distan t galaxies at redshifts z = 5.7 (|Hu et alj|2004h and 
2 = 6.5 (|Taniguchi et all 120051 ) . Long-slit spectral observa- 
tions of massive lensing clusters have been employed to in- 
vestigate the faint end of the Lya em ission luminosity fu nc- 
tion over the interval 4.5 < z < 5.7 |Santos et alj|2004h in 
addition to generatin g samples of candi date systems in the 
interval 7 < z < 10 (|Stark et al J 120071 ). In comparison to 
narrow-band observations, long slit spectroscopic observa- 
tions of Lya emitting sources typically probe a lower back- 
ground per resolution element yet are normally restricted to 
relatively small volume studies due to the limited field of 
view of a spectrograph slit. 

A critical unknown factor determining the visibility of 
the Lya emission feature in z > 7 galaxies is the physical 
state of the intervening IGM. Absorption studies of z ~ 6 
quasars appear to hav e identified the v ery end of the global 
re-ionisation process l|Fan et al.l 120061 ) while measurement 
of the optical depth of electron scattering at large angular 
scales in the CMB points to a typical (t hough not defin i- 
tive) epoch of reionisation around z — 11 (|Page et alj |2006). 
Clearly then, galaxies at redshifts 7 < z < 11 may be lo- 
cated in a partly ionised IGM where the local fraction of 
neutral hydrogen around individual galaxies may be suffi- 
cient to attenuate the Lya emission signature. However, in- 
dividual Lya emitting galaxies have been identified in deep, 
narrow z-band sur veys and confir med spectrosco pically at 
redsh ifts z = 6.56 l|Hu et aljfeool ) and z = 6.96 (|lve et all 
2006). Prompted by such observations, numerous theoret- 
ical studies have been undertaken to compute the escape 
fraction of ultra-violet (UV) photons from a volume of high- 
redshift intergalactic hydrogen ionised as a result of star 
formation occuring in an embe dded galaxy (lHaimanll2002i ; 
ISantosll2004l ; lBarton et al.ll2004h . While the detailed proper- 
ties of the Lya line transmitted through such a medium are 
necessarily model dependent, e.g. depending upon the mass, 
metallicity, star formation rate and initial mass function of 
the burst and the local density of the IGM, a range of plau- 
sible scenarios exist whereby a Hll region of sufficient size 
is created such that transmission of a partially attenuated 
Lya line occurs. 

The above considerations serve as the motivation for 
a NIR search for Lya emitting galaxies at z > 7. The 
search technique employs a narrow J-band filter centred 
on 1.187(i.m and is sensitive to the signature of Lya emit- 
ting galaxies located at a redshift z = 8.8 (termed z ~ 9 
in the following text). The remaining sections are organ- 
ised as follows: in Section [2] we describe in further detail the 



construction of the narrow-band survey. In Section [3] we de- 
scribe the techniques used to process the data and identify 
candidate z ~ 9 Lya emitting galaxies. Finally in Section 
[4] we determine the co- moving volume at z ~ 9 sampled in 
terms of the Lya emission luminosity and compare this to 
a reasonable range of expected properties of z ~ 9 galaxies. 
Throughout this paper, values of Om.o = 0.3, Qa,o = 0.7 and 
Ho — 70kms _1 Mpc _1 are adopted for the present epoch 
cosmological parameters describing the evolution of a model 
Friedmann-Robertson- Walker universe. All magnitude infor- 
mation is quoted using AB zero point values. 



2 THE ZEN SURVEY 

The desirability of detecting very high redshift Lya emit- 
ting galaxies forms the motivation for the "z equals nine" 
(ZEN) survey. Lya emitting galaxies occupying narrow red- 
shift intervals at z > 7 will present a characteristic emission 
excess signature in infrared photometry employing a com- 
bination of narrow and broad band filters. Narrow-band, 
NIR filters tuned to sample regions of night sky emission 
devoid of strong terrestrial line features provide access to 
relatively low background signals and thus permit sensitive 
imaging observations to be executed. In what we refer to as 
ZEN1 we constructed a 32 hour on-sky image of the Hubble 
Deep Field South (HDF-S) in the n arrow-band filter NB119 
employing the VLT/ISAAC facility (|Willis fc Courbinll2005l ; 
WC05 hereafter). Using deep, archival images of the field 
consisting of VLT/ISAAC J s -band and HST/WFPC2 opti- 
cal bands we were able to execute a sensitive test for faint 
narrow-band excess sources (i.e. NB < 25.2, J s — NB > 0.3) 
that remain undetected in optical bands - a practical defi- 
nition for candidate z ~ 9 Lya emitters. No candidate ZEN 
sources were identified in these observations. However, the 
study demonstrated that interloping low redshift emission 
excess sources could be successfully identified and rejected 
using deep optical images and that a detailed z ~ 9 vol- 
ume selection function could be computed in terms of Lya 
emission luminosity. 

As part of a new study - dubbed ZEN2 - we have ob- 
tained further narrow and broad band images directed to- 
ward the fields of three low redshift galaxy clusters: A1689, 
A1835 and AC114 (Table [T]). Each cluster acts as a gravi- 
tational lens and provides a spatially magnified view of the 
background universe. When considering unresolved sources 
the effect of this magnification is to increase the total bright- 
ness measured within a photometric detection aperture. We 
therefore use the presence of each cluster along the line of 
sight to provide a gravitational "boost" to the signal from 
putative z ~ 9 galaxies. The properties of each cluster have 
been described in detail in the literature and each possesses 
a well-determined gravitational lens model. Lens models de- 
scribing the clusters A1689, A1835 and AC114 are presented 
respectively in Limousin et al. (2007), Smith et al. (2005) 
and Campusano et al. (2001). Infrared observations of the 
three clusters were obtained with the VLT/ISAAC facility 
as part of ESO programmes 070.A-0643, 071.A-0428, 073.A- 
0475 and are summarised in Table [1] Optical observations 
of the three clusters consist of HST/ACS F850LP mosaics. 
These are described brie fly in Tablela nd in further d e tail in 
iBroadhurst et al l l|2005l ) for A1689 and lHempel et all (|2007f ) 
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Table 1. NIR data obtained for the three cluster fields. 



Cluster 


a (J2000) 


S (J2000) 


Redshift 




NB119 observations 






J-band observations 












■*V exp 


DIT (s) 


NDIT 


t-exp (s) 


AT exp 


DIT (s) 


NDIT 


texp (s) 


A1689 


13:11:30.1 


-01:20:17.0 


0.18 


54 


110 


3 


17820 


57 


35 


4 


7980 


A1835 


14:01:02.0 


+02:51:46.7 


0.25 


82 


100 


3 


24600 


10 


45 


3 


5400 


AC114 


22:58:47.7 


-34:48:04.1 


0.06 


78 


100 


3 


23400 


40 


45 


3 


5400 



Table 2. Properties of the HST/ACS F850LP images of each 
field. 



Cluster 




Image scale 
("/pix- 1 ) 


5cr limiting magnitude 
within 0'.'7 aperture 


A1689 


28600 


0.05 


27.48 


A1835 


18220 


0.04 


26.92 


AC114 


18368 


0.05 


26.98 



for A1835 and AC114. The NIR observations described in 
this paper all fall completely within the ACS mosaic areas. 



3 DATA REDUCTION AND ANALYSIS 

The NIR NB119 and J-band data were processed using tech- 
niques essentially identical to those described in WC05 and 
we summarize them briefly here. Imaging data were a) dark 
subtracted using standard calibration frames, b) corrected 
for varying pixel response using twilight sky exposures, c) 
sky-subtracted having masked array regions containing ob- 
jects above a specified ADU level, d) corrected for both high- 
and low-frequency spatial artefacts, and e) shifted to a com- 
mon pixel scale and coadded using a suitable pixel weighting 
and rejection algorithm. The image quality in each field and 
filter combination was computed as the mean full-width at 
half maximum (FHWM) of a sample of bright, stellar sources 
visible in each image and is displayed in Table 

The J-band observations were obtained under photo- 
metric conditions and were placed on an absolute flux scale 
using observations of standard stars taken with the science 
data. Reference magnitudes for standard stars observed with 
the NB119 filter are not available. The NB119 data were 
therefore placed on an absolute flux scale by comparing the 
flux measured within a 5" diameter circular aperture for a 
sample of bright, isolated stellar sources visible in both the 
J-band and NB119 images of each field. The colour term 
J — NB for these reference sources is approximately zero 
and the relationship between J-band reference magnitudes 
and NB119 instrumental magnitudes is linear and displays 
a gradient of unity. A total of 12, 16 and 16 such refer- 
ence sources were employed for the fields A1689, A1835 and 
AC 114 respectively. 

Source detection and photometry were performed on 
each image using the SExtractor software package (Bertin 
and Arnouts 1996). Once again, we continue the approach 
outlined in WC05, i.e. source detection is optimised for the 
detection of marginally resolved or unresolved objects and 
source fluxes are measured in circular apertures of diameter 
0'.'7. The required correction to convert aperture photometry 
to pseudo-total photometry (assumed to be a 5" diameter 
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Figure 1. Mean detection probability of simulated point sources 
as a function of AB magnitude within the NB119 (square sym- 
bols and dashed line) and J images (circular symbols and solid 
line). The horizontal dotted line indicates the 90% completeness 
threshold. 

circular aperture) was computed for bright stellar sources in 
each field. Corrections generated for each image are given in 
Table |3j though it should be noted that subsequent calcula- 
tions based upon image photometry all use the 0'.'7 aperture 
values. 

The magnitude corresponding to the 90% point source 
recovery threshold in each image (mgo) was computed by 
introducing artificial point sources into the reduced image 
of each field and determining the fraction recovered (see 
WC05) . This procedure produces a low resolution (3" pixels) 
image of the varying depth of each image. In addition, this 
process generates a map of photometric uncertainty across 
each field taking into account contamination by bright galax- 
ies. The average value of the completeness across each image 
is also displayed in Figure Q] as a function of source magni- 
tude. The magnitudes corresponding to the 90% complete- 
ness limit in each field and filter combination are displayed in 
Table The typical signal-to- noise ratio (SNR) of a source 
displaying mgo in each field is approximately 15. Note that 
for the purposes of computing the survey selection function 
in Section [3] we employ the two dimensional completeness 
information available for each field. 

Candidate z ~ 9 Lycv emitting galaxies are identified 
as narrow-band excess sources relative to the J-band refer- 
ence filter. Figure [2] displays the narrow-band excess versus 
NB119 magnitude for each of the target fields. We identify 
emission line sources as those displaying a positive J — NB 
signature in excess of the local 3a uncertainty in J — NB. At 
the 90% completeness limit of each field this corresponds to 
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Table 3. Image quality and photometric properties of each reduced image. 



Cluster 


Image quality (") 
NB119 J-band 


AB Depth at 90% completeness 

(0'.'7 aperture) 
NB119 J-band 


Correction to total magnitudes 
NB119 J-band 


Lya flux sensitivity 
(xl0~ 18 ergs s _1 cm~ 2 ) 
image plane source plane 


A1689 


0.47 


0.49 


24.5 


24.9 


0.63 


0.68 


3.7 


0.48 


A1835 


0.42 


0.42 


24.5 


24.8 


0.56 


0.55 


3.6 


2.38 


AC114 


0.45 


0.42 


24.5 


24.8 


0.63 


0.60 


3.8 


2.67 




NB119 (AB) 

Figure 2. Distribution of narrow-band excess J — NB versus NB119 magnitude for the three cluster fields. Values of J — NB versus 
NB119 for sources in each field are indicated by dots. For each panel the dotted curve indicates the three sigma uncertainty in the 
narrow band excess and the solid horizontal line indicates J — NB = O.The vertical dashed line in the each panel indicates the value 
NBII990 -Sources displaying a narrow-band excess greater than the three sigma uncertainty arc highlighted using a box. For clarity only 
sources satisfying NB > 22 or notable brighter sources are marked in this manner. Sources failing the above narrow-band excess measure, 
yet close to either the selection envelope or the magnitude cut, are indicated with a triangle. These sources were also inspected visually 
and excluded as potential ZEN sources due to their detection in the corresponding ACS z-band image. 



a colour excess J — NB ^ 0.3. Of the sources satisfying this 
narrow-band excess threshold, all are ultimately detected in 
HST F850LP images of each field. We therefore associate 
these sources with intervening emission line galaxies (e.g. 
[Oil] 3727, H/3 or Ha) located at redshifts that place the 
emission feature in the NB119 filter (c.f. WC05). The field 
of view of the NIR images of each cluster cover an area of 4 
square arcminutes in each case. Within this area the images 
of clusters A1689, A1835 and AC114 respectively contain 



4, 13 and 21 interloping emission line sources down to the 
observed frame magnitude limit mgo appropriate for each 
field. In addition, a small number of sources in each field 
are detected in the narrow-band but remain undetected in 
the accompanying J-band. All of these sources were investi- 
gated and ultimately associated with faint sources in optical 
HST observations. Therefore, no candidate z ~ 9 Lya emit- 
ting galaxies have been detected in any of the three fields 
studied. 
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4 PROBING THE LUMINOSITY FUNCTION 
OF z ~ 9 LAE GALAXIES 

The non-detection of candidate z ~ 9 Lya emitting galaxies 
in the three cluster fields may be understood by comput- 
ing the total volume sampled at each Lya luminosity. To 
achieve this the survey flux sensitivity as a function of solid 
angle must be corrected for the magnification introduced by 
the cluster lensing potential in each field, distance dimming 
and the effect of partial transmission of the Lya line by the 
NB119 filter. 

The gravitational potential associated with each galaxy 
cluster creates a magnified view of the background universe. 
This magnification can be considered as a spatially varying 
transformation between the source plane and image plane 
geometry of a particular field. The detection limits displayed 
in Figure[T]are associated with the image plane of each clus- 
ter. The corresponding source plane detection map f or each 
clust er was computed using the LENSTOOL package l|Kneibl 
1993) and a model potential describing each cluster. Each 
pixel in the detection map (location and area) was trans- 
formed to the source plane accounting for the varying lens 
deflection angle as a function of sky position before being re- 
assembled onto a uniform pixel grid. Images displaying the 
individual stages in this procedure are displayed in Figure 

m 

In Appendix A of WC05 we described how the mag- 
nitude limit of a particular narrow-band image could be 
transformed into a Lya emission line luminosity at z ~ 9 
by accounting for the equivalent width criterion applied to 
select narrow-band excess sources and the partial transmis- 
sion of the Lya line by the narrow-band filter. The mean 
Lya flux sensitivity toward each cluster field is given in Ta- 
ble [3] and is computed over both image plane and source 
plane (i.e. de-lensed) pixels. The volume sampled as a func- 
tion of Lya luminosity, V(LL y a)i is then computed as the 
integral over the differential co-moving volume element out 
to the maximum redshift at which a source displaying the 
specified luminosity would be detected. We apply the same 
procedure to determine V(Li,yc*) in the current study with 
only a minor modification: rather than compute the volume 
sampling based upon the average depth over each field, we 
compute the volume sampled per pixel in the source plane 
detection map describing each cluster. The volume sampled 
as a function of Lya luminosity for each cluster field is then 
computed as the contribution from individual detection map 
pixels, weighted by the solid angle of each pixel. The volume 
selection function for each cluster is displayed in Figure 2] 
Each curve is computed assuming a Lya emission line of 
rest frame velocity width a v = 50 kms (see WC05 for 
additional details). 

The inverse of the volume selection function for the ZEN 
survey is equal to the cumulative space density of z ~ 9 
Lya emitting galaxies sampled as a function of their emis- 
sion luminosity (Figure [5). The region above each curve in 
Figure [5] indicates the region of the cumulative luminosity 
function of putative z ~ 9 Lya emitting galaxies that can 
be ruled out as a result of the non-detection of bona-fide 
2 ~ 9 sources It is instructive to compare these limits to the 
cumulative luminosity function both of observed Lya emit- 
ting galaxies at redshift z = 6.6 (Kashikawa et al. 2006) and 
to z ~ 9 Lya emitting galaxies simulated within a semi- 
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Figure 4. Co-moving volume sampled brighter than a given Lya 
luminosity. The black curves indicate the co-moving volume sam- 
pled toward each of the three cluster fields in this paper: A1689 
(dashed), A1835 (dotted) and AC114 (solid). The blue curves 
indicate the total co-moving volume sampled toward the three 
ZEN2 cluster fields (dashed) and the co-moving volume sampled 
toward the ZEN1 field (HDF-South) described in WC05 (solid). 
The red curve indicates the total co-moving volume sampled by 
the ZEN1 and ZEN2 surveys to date. 



analytic model of galaxy formation (Le Delliou et al 2006). 
We consider the implications for each population in turn. If 
we use the observed population of Lya emitting galaxies at 
z — 6.6 as a model for emission at z ~ 9 then the current 
areal coverage of the ZEN survey would have to be increased 
by a factor of at least three in order to provide a realistic 
constraint on the putative z ~ 9 population. The prospect 
for extending deep, narrow-band surveys at z ~ 9 to wider 
areal coverage is promisin g given the adv ent of both the DA - 
ZLE (|Horton et al.ll2004) and HAWK-I l|Casali et alj|2006h 
NIR cameras at the ESO VLT. If instead the properties of 
Lya emitting galaxies at z ~ 9 are described by the semi- 
analytic model of Le Delliou et al. then the prospects for 
their detection is less certain. A small region of the model 
described by a UV photon escape fraction of 0.2 has already 
been tentatively ruled out by the ZEN survey. However, re- 
ducing the escape fraction to 0.02 results in a proportionate 
decrease in the Lya luminosity and the detection of such 
a population with either DAZLE or HAWK-I will remain 
challenging. An interesting alternative approach employs the 
current generation of relatively wide field NIR cameras op- 
crating on 4m class telescopes. In Figure [5] we display the 
anticipated results of ZEN3, a narrow-band search for Lya 
emitting galaxies at z ~ 8 currently underway using the 
Ca nada France Hawa ii Telescope (CFHT) WIRCam facil- 
ity (|Puget et al.ll2004T ). The exceptional volume sampling of 
such wide field cameras will permit a very sensitive test of 
the space density of putative z ~ 8 emitters, whether based 
upon observed z — 6.6 or model populations. 



6 J. P. Willis, F. Courbin, J. -P. Kneib, D. Minniti 




A1835 



• • 


■ 


• - 1 




J - 


■ ■ 1 


: ; j ./ *\ ■ > 
• • * i . . * W 
■ . WHS ^ 


■ 

mm ■ / 






■ -L 


»■ 




m 


■ i 


AC114 ■ < ■ • 


m 


^ ■ 



G 




i 



1 "J 




Figure 3. Grcyscale images demonstrating the processing applied to the two dimensional sensitivity maps. All images measure two 
arcminutes on a side with North up and East left. Panels for each cluster are arranged in rows. Left panel: narrow-band image of each 
field. Centre panel: image plane source detection sensitivity, mgo (see text for more details). Lighter regions indicate fainter sensitivity 
levels. Right panel: Lya flux detection sensitivity computed for a source plane located at z = 8.8. Lighter regions indicate fainter 
sensitivity levels. The field distortion arises from the inversion of the image plane sensitivity map though the cluster potential. In each 
panel the red and yellow contours indicate respectively the critical and caustic lines corresponding to a source located at z = 8.8. 
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Figure 5. A comparison of the limits imposed upon the z > 7 
Lya luminosity function by existing and planned NIR narrow- 
band surveys. The blue curves indicate the region of the cumula- 
tive space density versus Lya luminosity sampled by observations: 
the ZEN1+2 surveys (solid), the wide area, shallow depth ISAAC 
survey of Cuby et al. (2006) (dotted) and the planned sensitiv- 
ity of ZEN3, a wide area NB survey employing CFHT WIRCam 
(dashed). The point s indicate the ob s ervat ions of z = 6.6 Lya 
emitting galaxies by Kashikaw a et al ] j2006t) : circles indicate the 
photometric sample corrected for completeness while the trian- 
gles indicate the spectroscopically confirmed sample. The red 
curves indicate the Lya luminosity function simulated at z = 9 by 
iLe Delliou et al. assuming a Lya escape fraction of 0.02 

(solid) and 0.2 (dashed) respectively. 
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